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I.  INTRODUCTION 


A quantitative  description  of  the  interaction  of  electrons  with 

solids  over  a wide  range  of  energies  is  a subject  of  importance  in  a 

wide  variety  of  basic  and  applied  physical  problems.  Theoretical 

calculations  of  energy  loss  and  range  of  electrons  in  many  materials 
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have  formed  the  basis  of  at  least  two  extensive  tabulations.  ’ Both 
of  those  tabulations  are  based  on  the  Bethe  theory  of  stopping  power 
and  are  restricted  to  electron  energies  =^10  keV.  To  complement  these 
results  we  have  employed  several  theoretical  models  to  provide  calcu- 
lations of  inverse  mean  free  path,  energy  loss,  csda  range,  and  straggling 
for  electrons  with  energies  < 10  keV.  Tables  of  these  quantities  are 
now  available  for  the  solids  A1  and  A^Og  (Reference  3);  Si  and  Si 0^ 

(Reference  4);  Ni , Cu,  Ag,  and  Au (Reference  5);  and  Ge  and  GaAs  (Reference  6). 
These  tables  should  provide  useful  guides  for  interpretation  of  experi- 
mental data  as  well  as  input  for  calculations  in  applied  areas. 

The  work  presented  here  for  the  organic  insulator  polystyrene, 

o 7 

-(CgHg)-,  employs  a model  insulator  theory  ’ ’ to  describe  the  response 
of  the  valence  band  electrons.  The  states  of  the  tightly  bound  K-shell 
electrons  are  assumed  to  retain  a free-atom-like  character  so  the 
excitation  of  these  electrons  to  the  continuum  is  described  by  cross  sections 

O 

derived  from  atomic,  generalized  oscillator  strengths  (GOS's).  In  the 
following  sections  we  describe  the  calculation  of  differential  inverse 
mean  free  paths  (DIMFP's)  for  interaction  of  an  electron  with  the  valence 
band  or  carbon  K-shell  electrons  in  the  solid  polystyrene  and  the  derivation 
of  inverse  mean  free  path  (IMFP)  and  energy  loss  from  these  DIMFP's. 


4 


Results  are  presented  graphically  for  the  IMFP  and  energy  loss  (or 
stopping  power  of  the  polystyrene)  and  in  tabular  form  for  the  IMFP, 
stopping  power,  csda  range,  and  straggling  for  electrons  of  energy  from 
a few  electron  volts  through  10  keV. 

II.  GENERAL  FORMULATIONS 

A charged  particle  passing  through  a solid  interacts  with  a large 
number  of  electrons  simultaneously,  and  it  is  thus  appropriate  to  speak 
of  a mean  free  path  of  the  charged  particle  for  energy  transfer  to  the 
solid.  Assuming  the  effect  of  the  charged  particle  on  the  medium  may 
be  described  in  first  Born  approximation,  the  inverse  mean  free  path, 
differential  in  momentum  transfer  Ilk,  and  energy  transfer  liu),  for  a 
particle  of  velocity  v is  given  by 
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where  e(k,u>)  is  the  dielectric  response  function  of  the  solid.  * We 
assume  in  this  work  that  the  solid  is  Isotropic  and  homogeneous  so  that 
c is  a scalar  function  of  k and  u>. 


For  our  calculations  of  Inverse  mean  free  path,  stopping  power, 

etc..  It  Is  sufficient  to  compute  inverse  mean  free  paths  differential 

In  energy  transfer  only.  This  differential  Inverse  mean  free  path 

2 

(DIMFP)  for  energy  loss 'flu)  by  an  electron- with  energy  E ■ mv  /2  In  the 
solid  Is  given  by 


x(E,flu>)  = 


dffcy = + ir  Im[r[ra] 


— , 2 2 
where  flk+  ; /5m  [.f  + /E-Aw]  and  aQ  = A /me  . This  expression  assumes 

that  the  energy -momentum  relation  for  a swift  electron  in  the  solid  does 

not  differ  appreciably  from  that  of  a free  electron  in  vacuum. 

Given  e(k,w)  for  the  solid,  the  quantities  of  interest  here  follow 

directly  from  x(E,Au>).  The  inverse  mean  free  path  of  the  electron,  p, 

is  given  by  integrating  over  allowed  energy  transfers  as 


v.(E)  = /dOv)  t(E,1W). 


The  rate  of  energy  loss  of  the  electron,  or  the  stopping  power  of  the 
medium.  Is  given  by 


S(E)  3 -dE/dx  *Jd(1ta)  flu.  x (E  ,flu>) , 


and  the  mean  square  energy  loss  per  unit  path  length  by 


G2(E)  h J d (flu.) (flu.) 2 x(E.flu)). 
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With  these  results  we  may  calculate  the  range  of  an  electron  in  the 
continuous-slowing-down  approximation  (csda  range)  by 


Ro(E)  =JE  dEVS(E’) 


(6) 


and  the  mean  square  fluctuation  in  the  range  or  "range  straggling" 
will  be  calculated  from  Eqs.  (4)  and  (5)  as11 


(R  - R 0)aV  s J dE1  «2(E')/[S(E' )]3 


(7) 


For  our  tabulations  we  take  the  lower  limit  in  the  integrations  of 
Eqs.  (6)  and  (7)  as  Eq  = 10  eV. 


III.  DIMFP  FOR  THE  VALENCE  BAND 


The  model  insulator  theory  used  to  derive  the  DIMFP  for  inter- 
action of  an  electron  with  the  valence  band  electrons  has  been  described 

3 6 7 

and  employed  in  several  previous  calculations.  ’ Instead  of  repeat- 
ing the  detailed  formulae  here,  we  present  graphically  some  of  the  steps 
required  to  obtain  the  energy  loss  function,  Im(-l/e(k,w)J , which  is  the 
key  ingredient  in  the  calculation  of  the  DIMFP. 

The  first  step  in  applying  the  model  insulator  theory  is  to  fix 
the  adjustable  parameters  by  fitting  the  theoretical  expression  for  the 
imaginary  part  of  the  dielectric  response  function  in  the  optical  limit  (k*o) 
to  experimentally  determined  values  of  this  quantity.  In  Figure  1 we 
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show  the  results  of  the  fit  using  experimental  values  obtained  by  Inagaki  ; 

12  I 

et  al.  for  polystyrene.  The  valence  band  is  assumed  to  result  from  a 

combination  of  three  ground  state  orbitals  and  the  fit  shown  in  Figure  1 

gives  the  following  set  of  parameters  (as  defined  in  References  3,  6,  and 

7): 


i 

^0 

■nu,Bi(ev) 

"1 

1 

0.29 

5.80 

2.5 

2 

0.87 

8.50 

31.5 

3 

1.90 

25.0 

7.5 

where  n^  is  the  number  of  valence  electrons  per  monomeric  unit  accounted 

for  by  each  level  and  6 = 1/2  for  each  of  the  three  levels.  Note  that 

we  account  for  41.5  electrons  per  monomeric  unit  in  the  valence  band  instead 

of  the  expected  40.  This  redistribution  of  electron  numbers  between  the 

core  and  valence  electrons  is  due  to  oscillator  strength  couplinq  between  the 

1 2 

core  and  valence  levels.  Since  there  are  56  electrons  per  monomeric 
unit,  we  will  account  for  an  effective  number  of  14.5  carbon  K-shell  electrons 
(or  7.25  K shells)  per  monomeric  unit  in  our  calculations  of  DIMFP's 
in  Section  IV.  The  density  of  polystyrene  for  these  calculations  is  taken 

3 

to  be  1.05  g/cm  . With  the  molecular  weight  of  104.14  g/mole  for  polystyrene, 

-3  °3 

this  density  corresponds  to  6.07  x 10  monomeric  units/A  . 

As  a further  comparison  we  show  in  Figure  2 the  enerqy  loss  function  in 

12 

the  optical  limit  calculated  from  the  experimental  data  and  calculated  from 
the  model  insulator  theory  using  the  parameters  determined  above.  Quite 
reasonable  agreement  is  seen  in  both  Figures  1 and  2. 

The  extension  of  the  energy  loss  function  to  arbitrary  values  of 
momentum  transfer  as  determined  by  the  model  insulator  theory  is  illustrated 
in  Figure  3 where  energy  and  momentum  transfer  are  given  in  atomic  units. 
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Calculations  of  D 1 MFP  from  Eq.  (2)  are  illustrated  in  Figure  4 where 
we  plot  E : as  a function  of  energy  transfer  for  several  values  of  electron 
energy,  with  all  quantities  expressed  in  atomic  units. 

IV.  DIMFP  FOR  CARBON  K SHELL 


Frcm  a general  expression  for  the  dielectric  function  of  a 
homogeneous,  isotropic  system1''  we  may  show  that  for  values  of  u>  which 
correspond  to  ionization  of  a given  inner  shell  in  a solid  that 

2 

Im  M/e(k.u.)]  * Im  e(k,w)  % , (8) 

lliui  Qui 


where  df/dw  is  the  GOS  and  n is  the  number  of  those  inner  shells  per  unit 
volume  in  the  solid.  Equation  (2)  thus  leads  to 


i (E.tu) 


8,,ao2n  [K+  dk 

(E/R)  (1WR)  K k 


df(k,m) 

d(flw) 


(9) 


where  R = e^/2aQ  = 13.6  eV. 

Generalized  oscillator  strengths  for  ionization  of  electrons  from 

O 

the  K shell  of  carbon  have  been  calculated  by  McGuire.  These  GOS 
values  have  been  used  in  Eq.  (9)  to  obtain  the  differential  cross 
section  do/d(lL)  i/n.  Some  typical  results  are  shown  in  Figure  5 
for  several  values  of  electron  energy.  The  binding  energy  of  the  K-shell 
electrons  in  carbon  is  ^ 282  eV.  As  discussed  in  Section  III,  we  account 
for  an  effective  7.25  K shells  per  monomeric  unit  in  calculating  r. 
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V.  tXl'MANGl  CORRLCltn  DIMrr's  and  rORMULAf  I OB  THf  TARUI  AT  10NC- 


Wo  have  Included  the  effect  of  electron  exchange  in  our  calcu 
lations  in  a simple  manner  Rased  on  the  form  of  the  Mott  formula 
(nonrelat ivistlc  Miller  fonmila)  for  scattering  of  an  incident 
electron  with  a free  electron.  The  cross  section  for  finding  a 
scattered  electron  with  energy  W per  unit  energy  interval  is  given 


4 

dt*  ne 

aw  * nr 


i + jl  . ! 

W‘  (E-WV  W(l-W) 


(10) 


for  an  incident  electron  of  energy  I . except  for  energies  close  to 
W * 0 and  W * E.  Near  W = 0 and  W * l the  interference  term  (third 
term  on  the  right  side  of  Eq.  (10)) is  small  compared  with  the  first 
or  second  term,  respectively. 

The  D1MFP  for  excitation  of  an  electron  from  a particular  state 
i may  be  written  in  the  form 

^(E.Hm)  * | F.(t.1T„0.  (11) 

If  wo  assume  that  the  width  of  the  level  from  which  an  electron  i>- 
excited  is  quite  narrow,  we  obtain  from  lq.  (11)  the  IM  Ml  I'  for  production 
of  a secondary  electron  with  energy  I ^ as 

- ! Mt.r'Sr ) , (ir) 


B t fi 

where  E^  is  the  binding  energy  of  the  1 level  (a  positive  quantity) 
The  exchange  corrected  DIMfR  is  taken  as 


(E.liui)  = — I F ^ ( E (E  ,E+E^-fi(o) 


•/E  Fi(E.‘fta,)Fi(E.E+E®-1i(o)  *| 


• (13) 


Since  Ei.  « l/(fL>)‘  for  large  E and  fL,  Eg,  (13)  reduces  in  this  limit 
to  the  form  given  by  Eq.  (10).  The  factor  1 - reduces  the 

contribution  of  the  third  term  in  Eq.  (13)  as  E ♦ E.*?.  This  form  for 
the  exchange  corrected  DIMFP  has  been  used  in  our  calculations  for 
the  inner  shell  and  for  the  valence  bands  (since  our  model  assumes  the 
width  of  these  levels  to  be  quite  narrow). 

If  we  now  define  the  more  energetic  of  the  two  electrons  after 
collision  to  be  the  primary  and  account  for  exchange  through  Eg.  (13), 
Eq.  (3)  gives  the  contribution  to  the  inverse  mean  free  path  due  to 
excitation  of  an  electron  from  the  i^1’  level  as 


'i<E>  • f 

J n 


(E+E<?)/2 


d (ffu> ) ijW'(E,1i..i). 


Similarly,  for  the  stopping  power  and  mean  square  energy  loss  per  unit 
path  length,  we  have  from  Eg.  (4)  and  Eg.  (5) 


(E+El?)/2 


S i ( E ) 


d(fL)  ft-o  ,®XC  (E.-M 


lb 


and 


07(E)  = 


/ 


(E+E“)/2 


d (-M  ('Ku>)‘ 


exc 

[i 


(E.-ffu,). 


For  the  remaining  calculations  we  form  the  sums 
Sexc<E>  - J (E) 


06) 


(17) 


and 


‘ic>E> 


= z 


^(E), 


(18) 


where  the  index  i includes  the  terms  appropriate  for  a given  solid, 
including  exchange  corrections  as  indicated  above.  The  csda  range  is 
calculated  from 

Rom  */E  dE7sexc(E')  (19) 

10eV 

corresponding  to  an  electron  slowing  down  in  a continuous  manner  from 
an  energy  E to  10  eV.  The  mean  square  fluctuation  in  the  csda  range 
based  or.  Eq.  (7)  is  calculated  as 

(R-Ro>AV  -/E<iE'  "exc  <E')/tSexc(E')]3.  (20) 

10eV 
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VI.  RESULTS 

Before  presenting  the  tabulations  for  polystyrene,  we  discuss 
briefly  the  results  for  IMFP  and  stopping  power.  In  Figure  6 the  IMFP's 
are  shown  for  interactions  with  the  valence  band  (sum  of  contributions 
from  the  three  levels,  Section  III)  and  with  the  carbon  K shell.  The  K- 
shell  contribution  is  quite  small  compared  to  the  valence  band  contribution 
and  amounts  to  ^ 1%  of  the  total  IMFP  in  the  energy  range  covered  here. 

We  have  found  no  measurements  of  electron  mean  free  paths  in  polystyrene 

1 A 

for  E ^ 10  keV.  However,  Swanson  and  Powell  ^ performed  characteristic 
energy  loss  measurements  using  20  keV  electrons  on  polystyrene  films  and 
determined  mean  free  paths  for  the  7 eV  and  21  eV  losses.  These  losses 
correspond  to  the  peaks  seen  in  the  energy  loss  function  derived  from 
optical  data  shown  in  Figure  2.  They  determine  a mean  free  path  x,  in 
A,  of  17,400  t 5,500  for  the  7 eV  loss  and  410  * 80  for  the  21  eV  loss. 

If  we  assume  these  represent  the  dominant  inelastic  loss  processes,  then 

o 

the  total  mean  free  path  0/tjQjAL  ^x7eV  + 1/,V21eV^  is  400  ^ ^ 20*) 

at  20  keV.  Extrapolating  our  IMFP  values  to  20  keV  yields  y v.  2.5  x 10  ''  A ' 
or  x = 1/u  = 400  A in  excellent  agreement,  possibly  fortuitous,  with  the 
experimental  result. 

In  Figure  7 we  show  the  contributions  to  the  stopping  power  of 

polystyrene  for  electrons.  Also  shown  is  the  stopping  power  derived  from 
1 ? 

Bethe-Bloch  theory  ’ for  E >_  10  keV.  At  10  keV,  the  Bethe-Rloch 

0 3 

result  is  S = 0.237  eV/A  for  polystyrene  of  density  p = 1.05  g/cm  . 

O 

Our  value  for  S = S (valence)  + S (K-shell)  is  S = 0.238  eV/A  which 
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VIII.  TABLES 


V “ ™ 


The  results  of  the  calculations  for  polystyrene,  at  a density  of 
3 

1.05  g/ctn  , are  given  in  the  following  tables: 

Table  1 presents  total  inverse  mean  free  path  and  contributions 
to  IMFP  due  to  interactions  with  valence  band  electrons  or  carbon 
K-shell  electrons,  Eq.  (14),  in  units  of  A for  incident  electrons 
with  energies  5 eV  v [ < 10  keV. 

Table  2 presents  total  stopping  power  and  contributions  to  the 
stopping  power  due  to  interactions  with  valence  band  electrons 

O 

or  carbon  K-shell  electrons,  Eq.  (15),  in  units  of  eV/A  for 
electron  energies  5 eV  <_  E 10  keV. 


Table  3 presents  the  csda  range,  Eq.  (19),  in  units  of  A,  the 

mean  square  energy  loss,  Eq.  (16),  in  (eV)  /A,  the  mean  square 

°2 

range  fluctuation,  Eq.  (20),  in  A1",  and  the  relative  range 

2 -i1' 

straggling  given  by  [(R-RQ)^y]'  /R0  , for  electron  energies  15  eV 
< E < 10  keV. 
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